Deep Photometry of GRB 041006 Afterglow: Hypernova Bump 

at Redshift z = 0.716 1 



K. Z. Stanek 2 , P. M. Garnavich 3 , P. A. Nutzman 2 , J. D. Hartman 2 , A. Garg 2 , 
K. Adelberger 4 , P. Berlind 5 , A. Z. Bonanos 2 , M. L. Calkins 5 , P. Challis 2 , B. S. Gaudi 2 , 
M. J. Holman 2 , R. P. Kirshner 2 , B. A. McLeod 2 , D. Osip 6 , T. Pimenova 3 , T. H. Reiprich 7 , 
in"; W. Romanishin 8 , T. Spahr 2 , S. C. Tegler 9 , X. Zhao 3 

o : 

{Sj , kstanek@cfa.harvard.edu, pgarnavi@miranda.phys.nd.edu, pnutzman@cfa.harvard.edu, 

?h jhartman@cfa.harvard.edu, agarg@cfa.harvard.edu, kurt@ociw.edu, 

pberlind@cf a . harvard . edu , abonanos@cf a . harvard . edu , mcalkins@cf a . harvard . edu 
pchallis@cf a . harvard . edu , sgaudi@cf a . harvard . edu , mholman@cf a . harvard . edu , 
■ rkirshner@cfa.harvard.edu, bmcleod@cfa.harvard.edu, dosip@lco.cl, 

tpimenov@nd.edu, thomas@reiprich.net, wjr@nhn.ou.edu, tspahr@cfa.harvard.edu, 

CN| ■ 

Stephen.Tegler@nau.edu, xzhao@nd.edu 

o\ : 

m ■ ABSTRACT 

(N ■ 

o : 

We present deep optical photometry of the afterglow of gamma-ray burst 
(GRB) 041006 and its associated hypernova obtained over 65 days after detec- 
P H i tion (55 i?-band epochs on 10 different nights). Our early data (t < 4 days) joined 

O ' with published GCN data indicates a steepening decay, approaching F v oc t~ 6 

& ■ at early times (t <C 1 day) and F v oc t~ 13 at late times. The break at 

tf, = 0.16 ± 0.04 days is the earliest reported jet break among all GRB after- 
glows. During our first night, we obtained 39 exposures spanning 2.15 hours 
^ ■ from 0.62 to 0.71 days after the burst that reveal a smooth afterglow, with an 

o3 1 rms deviation of 0.024 mag from the local power-law fit, consistent with pho- 

tometric errors. After t ~ 4 days, the decay slows considerably, and the light 
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curve remains approximately flat at R ~ 24 mag for a month before decay- 
ing by another magnitude to reach R ~ 25 mag two months after the burst. 
This "bump" is well-fitted by a k-corrected light curve of SN1998bw, but only 
if stretched by xl.38 in time. In comparison with the other GRB-related SNe 
bumps, GRB 041006 stakes out new parameter space for GRB/SNe, with a very 
bright and significantly stretched late-time SN light curve. Within a small sam- 
ple of fairly well observed GRB/SN bumps, we see a hint of a possible correlation 
between their peak luminosity and their "stretch factor" , broadly similar to the 
well-studied Phillips relation for the type la supernovae. 

Subject headings: galaxies: distances and redshifts — gamma-rays: bursts - 
supernovae: general 

1. Introduction 

At least some gamma-ray bursts (GRBs) are produced by events with the spectra and 
light curves of core-collapse supernovae (SNe), as demonstrated decisively by two GRBs 
which occurred in 2003, GRB 030329/SN 2003dh (e.g. Stanek et al. 2003; Matheson et al. 
2003) and GRB 031203/SN 20031w (e.g. Malesani et al. 2004). At redshifts of z = 0.1685 and 
z = 0.1055, respectively, these are classical GRBs with the two lowest redshifts measured 
to date (although GRB 031203 was most likely sub-luminous compared to other classical 
GRBs). GRB 980425 at z = 0.008 was also associated with "hypernova" 1998bw (Galama et 
al. 1998), but the isotropic energy of that burst was 10~ 3 to 10~ 4 times weaker than classical 
cosmological GRBs which might place it in a unique class. 

In addition, late-time deviations from the power-law decline typically observed for opti- 
cal afterglows have been seen in a number of cases and these bumps in the light curves have 
been interpreted as evidence for supernovae (for a recent summary of their properties, see 
Zeh, Klose, & Hartmann 2004). Possibly the best case of a supernova bump was provided 
by GRB 011121 which occurred at a redshift of z — 0.36 and thus would have had a rela- 
tively bright supernova component. A bump in the light curve was observed both from the 
ground and with the HST (Garnavich et al. 2003; Bloom et al. 2002). The color changes in 
the light curve of GRB 011121 were also consistent with a supernova, designated SN 2001ke 
(Garnavich et al. 2003). 

Due to this mounting evidence, there have been claims that all supernovae that produce 
GRBs are type Ic hypernovae. The bias toward type Ic is partly due to the perceived difficulty 
in producing a jet that escapes from a star with a massive envelope. But this prejudice may 
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not be justified; after all, it was recently believed that due to the large baryon content of 
supernovae it was not possible for any supernova to be a GRB source. It is therefore prudent 
to still assume that the range of SN types that are responsible for GRBs, and their properties, 
is an unsolved observational problem. In several cases there is evidence that GRBs could 
indeed be associated with other type of supernovae, such as type Iln (Garnavich et al. 2003), 
or normal type Ib/c or fainter hypernovae (Price et al. 2003; Delia Valle et al. 2003; Fynbo 
et al. 2004). Obtaining magnitudes, colors and spectra of more GRB supernovae is clearly a 
top priority in understanding the origin of long/soft bursts. 

GRB 041006 was detected and localized by the French Gamma Ray Telescope and the 
Wide Field X-Ray Monitor instruments aboard the High Energy Transient Explorer II at 
12:18:08 (UT is used throughout this paper) on 2004 October 6 (Galassi et al. 2004). It 
was classified as a "X-ray rich GRB", and it was similar to GRB 030329 in its lightcurve 
shape and spectral characteristics. Da Costa, Noel, & Price (2004) and Price, Da Costa, & 
Noel (2004) reported discovery of a new, fading optical source within the 5.0' HETE error 
circle, located at a = 00 h 54 m 50!2, S = +01 o 14'07" (J2000.0), and identified it as the GRB 
optical afterglow. Fugazza et al. (2004) and Price et al. (2004) obtained an absorption line 
redshift of z = 0.716 for the afterglow, which is fairly low for a GRB. Motivated by the con- 
tinued need to study the GRB/SN connection as discussed above, we undertook photometric 
monitoring of the optical afterglow of GRB 041006. In this paper, we present photometry 
of the afterglow of GRB 041006 obtained throughout the two months immediately following 
the burst, which resulted in a clear detection of a light curve bump, most likely due to the 
underlying hypernova. We discuss the properties of that hypernova compared to the sample 
of other well studied GRB-related SNe. 



2. The Photometric Data 

The photometric data we obtained are listed in Table l 1 . The majority of our data 
were obtained with the MegaCam CCD mosaic (McLeod et al. 2000) mounted on the MMT 
6.5m telescope. Additional data were obtained using the 1.8m Vatican Advanced Technology 
Telescope, the Magellan 6.5m Baade and Clay telescopes, and the Keck II 10m telescope 2 . 

All the data were reduced using DAOPHOT II (Stetson 1987, 1992). For the MMT 



1 The analysis presented here supersedes our GCN Circulars by Garnavich et al. (2004b) and Garg et al. 
(2004). 

2 All photometry presented in this paper are available through anonymous ftp on cf a-f tp . harvard . edu, 
in the directory pub/kstanek/GRB041006. Images are available by request. 
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GRB 041006 42" x 42" Clay 6.5m + MagIC 




UT 2004-1 1 -03 R-band 1 200 sec 



Fig. 1. — i?-band image of the afterglow [circle) of GRB 041006 obtained with the Clay 6.5m 
telescope on UT Nov. 3, 2004 (27.5 days after the burst). The brightness of the afterglow is 
R = 24.08. The photometric comparison star is marked with an "X" . This figure will not be 
included in the ApJL version of the published paper due to space constraints. 
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data, we have used 10-50 stars to establish the zero-point offset between frames. We found 
that a nearby star, marked "X" in Figure 1, is constant over many days to within 0.02 mag, 
with R = 21.90 using the calibration of Garnavich et al. (2004b) and Henden (2004). For 
other data, we use star "X" as the photometric calibrator. For the MegaCam data, where a 
Sloan r' filter was used, we found that it translates well into standard Rc filter magnitudes 
without a need for a color term. 



The transition between the phase dominated by the afterglow and the phase dominated 
by the supernova is marked by an abrupt change in the behavior of the light curve. In §3.1 
we consider photometry before this observed change, which occurs about 4.5 days after the 
burst. In §3.2 we consider data taken during the following two months. 



Figure 2 presents the i?-band light curve for GRB 041006. We have added some data 
published via the GCN (see the caption of Figure 2 for the references) to fill in the gaps in our 
data. All of the GCN data were brought to the common zero-point, and if no photometric 
error was reported, we adopted a conservative value of 0.2 mag. We have fit the GRB 041006 
data with the broken power-law model of Beuermann et al. (1999) (see also Stanek et al. 



where tf, is the time of the break, F^ is the i?-band flux at tb and s is a parameter which 
determines the sharpness of the break, where a larger s gives a sharper break. This formula 
smoothly connects the early time t~ ai decay (t <C tb) with the later time t~ a2 decay (t 3> 
U). Equation 1 has been used to describe the afterglow decay of e.g. GRBs 990510 (with 
s = 1) and 010222 (Stanek et al. 1999, 2001). The fit results in the following values for the 
parameters: «i = 0.59±0.13, a<i = 1.32±0.02, tb = 0.16±0.04 days (we have fixed the value 
of s = 2.5). Given the non-uniform data used (i.e. GCN plus our data), the quoted values 
for the errors of these parameters should be treated with caution. 

The resulting fit is shown as the dashed curve in Figure 2. The curve is a good fit 
to the early afterglow, with x 2 /DOF = 1.1. The fit yields a very early break time, tb = 
0.16 ±0.04 days, which is the earliest reported afterglow break (see the compilation of break 



3. Temporal Behavior in the Photometry 



3.1. Early Photometry: Days 1-4 



2001): 
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Fig. 2. — i?-band light curve for GRB 041006. Early data (days 1-4) are fit to a broken 
power-law (dashed curve). The fit yields a very early break time, tb = 0.14 days, which 
might be the earliest jet break reported to date. Filled points are our data, while open 
circles are data from GCN (Kinugasa & Torii 2004; Kahharov et al. 2004; Fugazza et al. 
2004; D'Avanzo et al. 2004; Monfardini et al. 2004; Misra & Pandey 2004; Covino et al. 
2004; Balman et al. 2004; Bikmaev et al. 2004). 
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times, tj e t in Table 1 of Friedman & Bloom 2005; see also Ghirlanda et al. 2004). Given the 
implications for the jet opening angle and the energetics of the burst, it would be valuable 
to further constrain the break time with robust calibration of other early data. 

Short-timescale variability has been seen in the light curves of several afterglows, in- 
cluding GRB 011211 (e.g. Holland et al. 2002) and GRB 021004 (e.g. Bersier et al. 2003). 
Motivated to look for such short-timescale variability in GRB 041006, the first night on the 
MMT we obtained 39 high S/N exposures of the afterglow, spanning 2.15 hours from 0.62 
to 0.71 days after the burst. No short-timescale variability is present in our data, and the 
rms deviation of 0.024 mag from the local power-law fit is consistent with random errors of 
the photometry. The afterglow of GRB 041006 joins other well observed afterglows, such as 
GRB 990510 (Stanek et al. 1999) and GRB 020813 (Laursen & Stanek 2003) in the category 
of very smooth afterglows. 



3.2. Later Photometry: Days 5-65 

Over the subsequent two months we obtained i?-band observations of the OT on a 
number of epochs. We combine our data with three additional i?-band photometric mea- 
surements published via GCN (Covino et al. 2004; Balman et al. 2004; Bikmaev et al. 2004). 
These observations are presented in Figure 3. 

As is apparent in Figure 3, the decay abruptly slows ~ 5 days after the burst. For 
slightly more than a month, the brightness then remains roughly flat at R ~ 24 mags. Over 
the next 30 days, the OT decays by another magnitude, falling to R ~ 25 mag in our final 
observations. 

This clear detection of a light curve "bump" strongly suggests the late-time dominance 
of the OT by an underlying supernova component. The supernova component peaks at 
roughly ~ 24 mag about 35 days after the burst, which is several magnitudes brighter than 
the extrapolated optical afterglow of the GRB at the time of this peak. 

To estimate the properties of the hypernova bump in GRB 041006, we correct for the 
effects of Galactic extinction using the reddening map of Schlegel, Finkbeiner, & Davis 
(1998). At the Galactic coordinates of GRB 041006, I = 124?74, b = -61?66, the foreground 
reddening is E(B — V) = 0.023, yielding an expected extinction of Ar = 0.06. We find that 
the OT light curve can be well modeled by a stretched SN 1998bw-like bump, k-corrected 
to redshift z = 0.716, combined with the continued power- law decay of the GRB optical 
afterglow. It is clear from the data that the amount of stretching of the SN 1998bw light 
curve needed is significant, and it is basically determined by the rapid decline by about a 
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GRB 041006 




Days after October 6.5126 

Fig. 3. — Late i?-band light curve for GRB 041006. Open circles are data from GCN, 
while the filled points show our data. The early-time broken power-law decay of the optical 
afterglow is shown as the dashed line. A k-corrected SN 1998bw light curve extended by 
a factor of 1.38 is shown as the dotted curve. The combined power law and stretched SN 
1998bw light curve is given by the solid curve. 
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magnitude from 36.9 to 56.5 days after the burst. The composition is shown as the solid 
line in Figure 3. Indeed, we find a good fit using as a template light curve of SN 1998bw, 
uniformly stretched in time by a factor of 1.38±0.06 (in addition to the 1+z cosmological time 
dilation), and with significantly brighter peak absolute magnitude (0.3 mag brighter). We 
did not fit for the brightness of the host, as its influence on the derived i?-band magnitudes 
would vary for our data set obtained with different telescopes, instruments and varying seeing 
conditions. To account for that, we adopt an asymmetric error bar for the peak brightness 
of the supernova of —0.1, +0.2 mag. 

Thus the hypernova component in the late-time OT of GRB 041006 is long-lasting, in 
contrast to e.g. GRB 01121/SN 2001ke (Garnavich et al. 2003) and GRB 030329/SN 2003dh 
(Lipkin et al. 2004; Deng et al. 2005), in which the bumps decay faster than SN 1998bw (see 
next Section for more discussion). As in the case of GRB 030329/SN 2003dh (Matheson et 
al. 2003), we find no need to introduce a time delay (either positive or negative) between 
the afterglow and the supernova component to fit the light curve data. In fact, adding 
a time delay parameter degrades the ^/DOF value slightly. The best fit time delay is 
—0.9 ± 2.7 days, i.e. the GRB and the SN were most likely simultaneous. 



4. Summary and Discussion 

We have obtained deep photometry of the optical transient associated with GRB 041006 
for the first two months following the burst. We find that the early-time i?-band light curve 
is well fit by a broken power law, with a pre-break index of a,\ = 0.59 and a steeper post- 
break index of «2 = 1.32. Our fit yields a break time of % = 0.16 days after the burst, the 
earliest jet break observed so far. Early observations also show a very smooth afterglow, 
with short-timescale deviations from the power-law decay on the same order (0.025 mag) as 
random errors in photometry. 

Continuing observations of the optical transient over the two months following the burst 
show an abrupt change from the t -1 ' 3 decay of the afterglow, consistent with the domination 
of the photometry by a SN component. The SN light curve peaks at a similar brightness to 
SN 1998bw but takes more time following the burst to evolve, by a factor of ~ 1.38. 

To put the light curve features of GRB 041006 into the context of previous GRB/SN ob- 
servations, we have adopted Figure 3 from Zeh, Klose, & Hartmann 2004 (hereafter ZKH04), 
which presents the deduced time stretch factors relative to SN 1998bw versus their de- 
duced luminosity ratios of GRB/SNe relative to SN 1998bw. They consider nine GRBs, 
observed before the end of 2002, for which there exists some evidence of a SN bump. They 
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Fig. 4. — Time stretch factor relative to SN 1998bw versus luminosity relative to SN 1998bw 
for six GRB/SNe (see also Figure 3 from ZKH04). Dashed lines at stretch factor = 1 and 
luminosity ratio = 1 represent SN 1998bw values. The four filled circles represent data for 
four GRB/SNe given in ZKH04. GRB 030329/SN 2003dh (triangle) and the SN component 
for GRB 041006 (star) are also shown (see text for more discussion). Plotted with the 
continuous line is the analogous relation for type la supernovae, adopted from Garnavich et 
al. (2004a). 
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have suggested a subdivision of these GRBs into a group of five with weak evidence for a 
bump, and a group of four with statistically significant evidence for a bump (GRBs 990712, 
991208, 011121, and 020405; see §3 of ZKH04). In our Figure 4, we present only this lat- 
ter group, with the additions of GRB 030329/SN 2003dh, which occurred after their cutoff 
date, and GRB 041006. For GRB 030329/SN 2003dh, we adopt the values for the parame- 
ters from a recent paper by Zeh, Klose, & Hartmann 2005 (hereafter ZKH05). It should be 
stressed that given the bumpy and long-lasting afterglow of GRB 030329, the exact param- 
eters for the underlying SN 2003dh are still subject to debate (Matheson et al. 2003; Lipkin 
et al. 2004; Deng et al. 2005). We decided not to include another recent GRB/SN case, 
GRB 031203/SN 20031w, as it was heavily obscured by foreground Galactic dust and resided 
in a rather luminous host galaxy. For that supernova, ZKH05 give values for the stretch of 
1.14 ± 0.16, and luminosity ratio of 1.65 ± 0.41. We note that, among this sample of six, 
GRB 041006 possesses the highest redshift, with z = 0.717. Also very interestingly, all GRBs 
with optical afterglows and measured redshifts less than that of GRB 041006 show evidence 
of late-time SN components (ZKH04). As apparent in Figure 4, GRB 041006 stakes out new 
parameter space for GRB/SNe, with a very bright and significantly stretched late-time SN 
light curve. 

The distribution of light curve stretch factors and peak brightnesses for hypernovae 
shown in Figure 4 is reminiscent of the relation found by Phillips (1993) for type la SNe 
(technique introduced by Perlmutter et al. 1997 for their sample of type la). We overplot the 
type la relation extended to fast declining events by Garnavich et al. (2004a) and extrapo- 
lated to large stretch factors. While the uncertainties for the hypernovae make a conclusive 
statement difficult, we note that both the hypernovae and type la events are likely powered 
by 56Ni decay and have similar peak luminosities implying comparable 56Ni masses. But 
with the present data available for them, hypernovae are still a long way from being useful 
distance indicators. More events like GRB 041006 will be possible to study soon, as the Swift 
satellite (Gehrels et al. 2004) is already providing accurate localizations for a good number 
of GRBs. 
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TV TTV Tm 

MMT 





.6657 


21 


.051 





.032 


120 


TV TTV TTT 

MMT 





.6675 


20 


.991 





.021 


120 


TV /TTV /TIT! 

MMT 





6693 


21 


,023 





.021 


120 


TV TH Tm 

MMT 





,6712 


21 


,070 





.023 


120 


TV TTV. Tm 

MMT 





,6730 


21 


,020 





,026 


120 


MMT 





,6749 


21 


,037 





,026 


120 


MMT 





,6767 


21 


,057 





,025 


120 


MMT 





,6785 


21 


,028 





.019 


120 


MMT 





,6804 


21 


,028 





,028 


120 


MMT 





,6822 


21 


,073 





,029 


120 


MMT 





,6840 


21 


,056 





,030 


120 


MMT 





,6859 


21 


,051 





,030 


120 


MMT 





.6877 


21 


.036 





.028 


120 


MMT 





.6895 


21 


.052 





.025 


120 


MMT 





.6913 


21 


.034 





.026 


120 


MMT 





.6931 


21 


.077 





.025 


120 


MMT 





.7014 


21 


.079 





.024 


120 


MMT 
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Table 1 — Continued 



A 1 


tic 


or 


I'exp 


Observatory 


0.7032 


21.050 


0.019 


120 


MMT 


0.7051 


21.044 


0.023 


120 


MMT 


0.7069 


21.047 


0.021 


120 


MMT 


0.7087 


21.052 


0.020 


120 


MMT 


0.7105 


21.036 


0.022 


120 


MMT 


0.7124 


21.086 


0.025 


120 


MMT 


1.7407 


22.330 


0.04 


1500 


VATT 


1.7851 


22.351 


0.020 


600 


MMT 


1.7931 


22.367 


0.022 


600 


MMT 


2.7341 


22.970 


0.05 


4800 


VATT 


2.7557 


22.878 


0.06 


600 


MMT 


2.9086 


22.980 


0.08 


600 


MMT 


2.9176 


23.072 


0.07 


600 


MMT 


3.7266 


23.330 


0.08 


4200 


VATT 


3.8531 


23.539 


0.07 


600 


MMT 


3.8613 


23.568 


0.07 


600 


MMT 


9.8622 


24.119 


0.16 


600 


MMT 


27.5763 


24.080 


0.10 


1200 


Clay 


30.7511 


24.040 


0.10 


1200 


Clay 


36.9040 


23.870 


0.11 


720 


Keck II 


56.5242 


24.770 


0.12 


1800 


Baade 


64.5516 


25.110 


0.23 


1200 


Baade 



a Days after 2004 October 6.5126 UT. 

b MMT: MMT telescope/MegaCam; VATT: Vati- 
can Advanced Technology Telescope; Clay: Magellan 
Clay teL/MagIC; Keck II: Keck II tel./ESI; Baade: 
Magellan Baade tel./IMACS 



